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Abstract 

The measurements of the Collins and Sivers asymmetries of identified hadrons pro- 
duced in deep- inelastic scattering of 160 GeV/c muons on a transversely polarised 
6 LiD target at COMPASS are presented. The results for charged pions and charged 
and neutral kaons correspond to all data available, which were collected from 2002 
to 2004. For all final state particles both, the Collins and Sivers asymmetries turn 
out to be small, compatible with zero within the statistical errors, in line with 
the previously published results for not identified charged hadrons, and with the 
expected cancellation between the u- and d-quark contributions. 
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1 Introduction 

To fully specify the quark structure of the nucleon at the twist-two level, the trans- 
verse spin distributions Axq(x) must be added to the better known spin-average distri- 
butions q(x) and to the helicity distributions Aq(x) [1]. Here x is the Bjorken variable, 
which represents the momentum fraction of the quarks inside the nucleon. The interpre- 
tation of the transversity distribution is similar to that of the helicity distribution, i.e. 
in a transversely polarised nucleon A?q is the difference of the number density of quarks 
with momentum fraction x and spin parallel or antiparallel to the nucleon spin. For a 
discussion on the notation see Refs. [2] and [3]. 

The distributions Axq are difficult to measure, since they are chirally odd and there- 
fore absent in inclusive deep inelastic scattering (DIS). They may instead be extracted 
from measurements of the single-spin azimuthal asymmetries in cross-sections for semi- 
inclusive DIS (SIDIS) of leptons off transversely polarised nucleons, in which a hadron is 
also detected in the final state. In these processes the measurable asymmetry is due to the 
combined effect of Axq and a chirally-odd fragmentation function (FF) which describes 
the spin- dependent part of the hadronization of a transversely polarised quark. At lead- 
ing twist, the existence of such a naively T-odd FF arising from final state interaction 
effects, was predicted by Collins [4] and is now generally known as the Collins effect. In 
the fragmentation of transversely polarised quarks it is responsible for a left-right asym- 
metry which is due to a correlation between the spin of the fragmenting quark and the 
transverse momentum p± of the produced hadron with respect to the quark direction. The 
Ri-dependent fragmentation function of a transversely polarised quark q into a spinless 
hadron h is thus expected to be of the form 

D T h q (z,p ± ) = D h q (z,p ± ) + A° T D h q (z,p ± ) ■ sin^ , (1) 

where D q is the unpolarised FF and the "Collins function" A^D q is the T-odd part of 
the FF, responsible for the left-right asymmetry, and z is the fraction of available energy 
carried by the hadron. Here ip is the difference of the azimuthal angles of the hadron 
transverse momentum and the quark spin, relative to the quark direction. As a result, 
in SIDIS off transversely polarised nucleons the Collins mechanism is responsible for a 
modulation in the azimuthal distribution of the produced hadrons given by 

N($ c ) = iVo(l + ec-sin$ c ). (2) 

In a gamma-nucleon reference system (GNS), in which the z-axis coincides with the virtual 
photon direction and the x-z plane is the lepton scattering plane, the "Collins angle" $c 
is $c = 4>h + <Ps — 7T. Here <ph is the azimuthal angle of the transverse momentum p£ 
of the outgoing hadron and cf>s is the azimuthal angle of the transverse spin vector S of 
the target nucleon, as shown in Fig. [TJ The measurable asymmetry ec is related to the 
convolution of the transversity parton distribution function (PDF) Aj-q and the Collins 
function. 

An entirely different mechanism was suggested by Sivers [5] as a possible cause 
of the transverse spin effects observed in pp scattering. This mechanism could also be 
responsible for a spin asymmetry in the cross-section of SIDIS of leptons off transversely 
polarised nucleons. Sivers' conjecture was based on a possible existence of a correlation 
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Figure 1: Definition of the azimuthal angle <ph of the transverse momentum p£ of the 
outgoing hadron and of the azimuthal angle 05 of the transverse spin vector S of the 
target nucleon. 

between the intrinsic transverse momentum k± of a quark and the transverse polarisation 
vector of the nucleon S, i.e. that the quark distribution q(x) could be written as 

q T (x, k±) = q(x, k ± ) + \S\ ■ A^q(x, k±) ■ sin tp', (3) 

where tp' is the difference of the azimuthal angles of the transverse spin of the nucleon and 
of the quark transverse momentum, relative to the nucleon direction. In SIDIS off trans- 
versely polarised nucleons the Sivers mechanism results in a modulation in the azimuthal 
distribution of the produced hadrons 

N(<f> s ) = N (l + e s -wi$s), (4) 

where the "Sivers angle" $5 = <fth — 4>s is the relative azimuthal angle between the 
transverse momentum of the hadron p\ and the nucleon target spin in the GNS. In this 
case, the measurable asymmetry is related to the convolution of the Sivers PDF A^q 
and the unpolarised FF D^. 

Since the Collins and Sivers terms in the transverse spin asymmetry depend on the 
two independent angles $c and $5, measuring SIDIS on a transversely polarised target 
allows the Collins and the Sivers effects to be disentangled, and the two asymmetries 
can separately be extracted from the data. The same is true for the other azimuthal 
asymmetries which appear in the general expression of the SIDIS cross-section in the 
one photon exchange approximation [6]. Correlations between the different terms can be 
introduced by a non-constant acceptance of the apparatus. 

Collins and Sivers modulations have been shown experimentally to be non zero 
by the HERMES measurements of pion asymmetries in SIDIS on a proton target [7,8]. 
Independent information on the Collins function has been provided by the azimuthal cor- 
relations in e + e~ — > hadrons measured by the BELLE Collaboration [9]. COMPASS has 
already published results for the Collins and Sivers asymmetries with a deuteron target 
for non-identified hadrons [10,11], and the details of the experimental technique can be 
found there. Here we present the Collins and Sivers asymmetries on a transversely po- 
larised deuteron target for identified hadrons, i.e. charged pions, and charged and neutral 
kaons, which put more stringent constraints on a flavour separated analysis of these new 
transverse spin effects. 
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2 The COMPASS experiment and the SIDIS event selection 

The COMPASS experiment is set up at the M2 beam line at CERN, using both 
muon and hadron beams. For this measurement a longitudinally polarised /i + beam of 
160 GeV/c momentum was scattered off a solid 6 LiD polarised target consisting of two 
cylindrical cells along the beam direction. The two cells were polarised in opposite direc- 
tions, so that data were taken simultaneously on oppositely polarised targets to reduce 
the systematic errors. The target polarisation direction could be set either parallel or 
orthogonal to the beam direction. 

Particle tracking and identification are performed in a two-stage spectrometer, cov- 
ering a wide kinematical range. A Ring Imaging Cherenkov detector (RICH-1) and two 
hadron calorimeters provide particle identification. The RICH-1 [12] detector is a gas 
RICH with a 3 m long C4F10 radiator covering the whole spectrometer acceptance. Two 
spherical mirror surfaces reflect and focus the Cherenkov photons on two sets of detectors 
outside the acceptance region. The photon detection utilises multiwire proportional cham- 
bers (MWPC) with segmented Csl photocathodes which detect photons in the UV region. 
The trigger system comprises hodoscope counters and hadron calorimeters. Veto counters 
installed in front of the target are used to reject the beam halo. A detailed description of 
the spectrometer can be found in Ref. [13]. 

The COMPASS Collaboration has taken data with the 6 LiD target polarised trans- 
versely with respect to the incoming beam direction in 2002, 2003, and 2004. The transver- 
sity data were taken for about 20% of the available beam time. During data taking, par- 
ticular care was taken to ensure the stability of the apparatus. The target polarisation 
was typically reversed every 5 days to reduce systematic effects due to the different ac- 
ceptances of the two cells. To guarantee that the acceptance of the detector fulfills the 
stability requirements the distributions of several physical quantities (x, p£, ■ ■ ■ ) have 
been monitored. 

In the data analysis, the selection of events required a "primary vertex" , defined by 
the incoming and the scattered muon tracks, and at least one hadron outgoing from the 
primary vertex. In addition the projected beam track was required to cross both target 
cells. Clean hadron and muon selection was achieved using the hadron calorimeters and 
considering the amount of traversed material. To select DIS events, cuts on the photon 
virtuality Q 2 > 1 (GeV/c) 2 and the mass of the hadronic final state W > 5 GeV/c 2 were 
applied. The requirement of 0.1 < y < 0.9, where y is the fractional energy of the virtual 
photon, limits the error due to radiative corrections and avoids contamination from ir 
decay (upper cut), and warrants a good determination of y (lower cut). To safely recon- 
struct the hadron azimuthal angle <frh a minimum transverse momentum p\ of 0.1 GeV/c 
is required. Furthermore, to select the current fragmentation region, a lower limit for the 
relative energy of the hadron z is required. In the following, all hadrons with z > 0.2 
define the "all hadron" sample. The analysis is performed also for "leading hadrons" only, 
for which we expect an enhancement of the physics signal albeit lower statistics. The 
"leading hadron" is defined as the most energetic hadron coming from the primary vertex 
and having z > 0.25. Detailed information on the data analysis can be found in Ref. [11]. 
Additional cuts specific to the analysis of the identified charged and neutral hadrons are 
described in the next sections. 

3 Charged particle identification 

Charged hadrons are identified as ir and K using the RICH-1 detector. Several 
variables have been monitored to ensure the stability in time of the RICH-1 response. 
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Figure 2: The Cherenkov angle as measured with RICH-1 versus the particle momentum 
as measured by the spectrometer. 



The monitored variables are the hit multiplicities of the MWPCs, the mass hypothesis 
likelihood (described later in this section), and the number of hadrons identified as ir and 
K normalised to the number of reconstructed tracks. Blocks of data are rejected if any of 
the variables deviates from the mean values by more than three standard deviations. 

The identification procedure relies on a likelihood function constructed with the 
photons detected in RICH-1 and associated to the charged particle trajectory. The like- 
lihood function uses the photons of the signal and the theoretical expectation from the 
Frank and Tamm equation, taking into account possible signal losses due to dead zones 
in the detector. The description of the background photons, coming from other particles 
in the event and from the beam halo, is taken from the photon detectors occupancy in 
the data. 

Since the number of Cherenkov photons depends on the velocity of the particle, for 
a given momentum the signal yield is different for different mass hypotheses. Thus for 
each track the likelihood is computed for different mass and the background hypotheses. 
The particle identification is made by choosing the mass hypothesis corresponding to the 
highest value of the likelihood. To assure a clear distinction from the background, a cut 
on the ratio of the highest likelihood to the background likelihood is made. To insure a 
good separation with respect to different masses, a cut on the ratio of the likelihood to 
the second highest likelihood is also made. The cuts on these variables have been tuned 
on subsamples of the data considering events in which at least two oppositely charged 
hadrons from the primary vertex have been found. The cuts for K identification have 
been tuned on the meson peak in the invariant mass distribution of the two charged 
hadrons, maximising the product of the <fi signal and the signal-to-background ratio values 
in the peak. To tune the cuts for pion identification, we used the same approach, but using 
the p peak in the invariant mass distribution. In addition, all these cuts have been verified 
on Monte Carlo data. 
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Figure 3: Purities of the pion (left) and kaon (right) samples as a function of x. The 
open (closed) points are for positive (negative) pions, while open (closed) squares are 
for positive (negative) kaons. The shaded area on the right picture gives the evaluated 
systematic uncertainty for the kaons, while the same is negligible for pions. 

In Fig. [2] the dependence of the Cherenkov angle measured by the RICH-1 against 
the particle momentum clearly indicates the different regions in which hadrons can be 
identified. In particular, the Cherenkov thresholds are visible: they are about 2.6 GeV/c 
for 7r, 9 GeV/c for K, and 17 GeV/c for protons. To assure a minimum number of detected 
photons the hadron momenta have to be 0.5 GeV/c above threshold for pions and 1 GeV/c 
above threshold for kaons. The values chosen correspond to 4 emitted photons. The upper 
limit for the momentum has been set to 50 GeV/c for both ir and K, corresponding to 
1.5 u mass separation between these two mass hypotheses. 

The purity of the identified hadrons, defined as the fraction of pions inside the 
identified pion sample, and the fraction of kaons inside the identified kaon sample, has 
been evaluated from the data. The results for this analysis are shown in Fig. [3] for positive 
and negative pions (left) and kaons (right) as a function of x. While the purity for pions is 
very high and almost independent from x, for kaons the purity depends on the momentum 
and the polar angle of the hadron, and as a consequence it increases with x. This trend 
is less pronounced as a function of z and given the smaller correlation with the polar 
angle. On average, purities are higher than 95% for pions, around 70% for negative kaons 
and around 80% for positive kaons. 

The final statistics for the different data taking periods are summarised in Table [TJ 
The use of the RICH-1 information was not possible for the 2002 data, since the func- 
tionality of the detector during the transversity run in that year was not satisfactory. 



4 Kg identification 

The reconstruction of the Kg relies on the detection of the two decay pions. The large 
acceptance of the COMPASS spectrometer provides a good efficiency for the detection of 
the pion pair. In this part of the analysis the RICH-1 detector is not used to identify the 
pions. 

The signature of K s events is a Vq vertex, i.e. a vertex with no incoming but two 
outgoing charged particles, where the two detected particles have opposite charge. The 
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Table 1: Final statistics for the "all hadron" and the "leading hadrons" samples. 









"all hadrons" 


(■10- 


1 


"leading 


hadrons" (-10 


-6) 


Year 


Period 




71 


K + 


K- 






71 


K + 






2002 


1 










0.021 










0.014 


2002 


2 










0.014 










0.009 


2003 




1.71 


1.49 


0.31 


0.20 


0.077 


1.10 


0.93 


0.24 


0.14 


0.052 


2004 


1 


1.54 


1.33 


0.27 


0.18 


0.063 


0.98 


0.82 


0.21 


0.13 


0.043 


2004 


2 


2.03 


1.76 


0.36 


0.24 


0.083 


1.30 


1.09 


0.27 


0.17 


0.056 


Total 




5.28 


4.58 


0.94 


0.62 


0.258 


3.38 


2.84 


0.72 


0.44 


0.175 



sum of the two outgoing particle momenta must point to the primary vertex and the 
invariant mass of the two particle system, assuming pion mass for each of them, must 
agree with the K° s mass. In order to identify K s from the primary vertex, all Vq vertices 
downstream of the primary vertex have been considered. Moreover, the outgoing tracks 
are not allowed to be additionally associated with any primary vertex and have to satisfy 
the criteria of section 2. To test the association of the secondary vertex to the primary 
vertex, the angle between the reconstructed momentum of the hadron pair and the vector 
connecting the primary and the secondary vertex is calculated. A maximum angle of 
10 mrad is accepted. To ensure a proper distinction of the primary and secondary vertices 
a cut on the distance of the primary and the secondary vertex is chosen on the basis of 
the signal-to-background ratio for the K s signal. A distance of about 10 cm yields a good 
background suppression. Fig. H] shows the Armenteros plot of the hadron pair, where the 
transverse momentum px of one of the hadrons relative to the hadron momentum sum is 
plotted vs. the difference of the longitudinal momenta over their sum {j>li — PL2) / (pli + 
Pl-i)- In this plot the K s band is clearly seen, as well as the A and A bands. To reduce the 
background due to e + e~ pairs a lower cut on of 25 MeV is applied. To finally identify 
the K s candidates a cut on the invariant mass is applied. The reconstructed invariant 
mass is required to be within ±20 MeV of the known mass as shown in Fig. [5J Since 
the width of the fitted peak is o ps 6 MeV, the region of ±20 MeV covers more than 
99 % of the signal. Finally, the kaon transverse momentum with respect to the virtual 
photon direction is required to be larger than 0.1 GeV/c to assure a good resolution in the 
measured azimuthal angle. After all these cuts, the signal-to-background ratio is about 
15, constant over all the kinematic range. 

Very much as for the charged hadron case, both "leading Kg" and "all K^" samples 
have been considered. In Table [1] the final statistics used for the asymmetry evaluation is 
given for all the periods. 



5 Extraction of the asymmetries 

The Collins and Sivers asymmetries for "all" and "leading" pions and kaons have 
been evaluated separately in kinematical bins of x, p\ and z while the other two variables 
were integrated over. The number of events N- k in the upstream and downstream target 
cell (k = u, d) for the two polarisations (+, — ) in a given bin (j = C, S refers to Collins 
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Figure 4: Armenteros plot of the hadron pair. The Kg band can be clearly seen as well 
as the A and A bands. 

and Sivers) can be written as 

N% = Ftn k aaf th {^) ■ (1 ± e% sin $,■). (5) 

Here F is the integrated incident muon flux, n the number of target particles, a the spin 
averaged cross-section and dj the product of angular acceptance and efficiency of the 
spectrometer. The quantities e^ fc are given by 

6 C,k = f ' P T,k ' D NN ■ Acoll , Cs,k = f ' P T,k ' A Siv , (6) 

where Acoii and A$i V are the "Collins" and "Sivers" asymmetries, related to the transver- 
sity and Sivers PDF respectively. The quantity D^n = (1 — £/)/(! — V + y 2 /2) is the 
transverse spin transfer coefficient from the target quark to the struck quark, / is the 
dilution factor and P^ k the absolute value of the polarisation of the target cells. The dilu- 
tion factor has been evaluated taking into account the radiative corrections for hadronic 
events, and is taken as constant, / = 0.38, known to 5%. The target polarisation P^ k has 
been measured for each cell and for each period [11] and averages to 48% with a relative 
error of 5%. 

From the measured number of events iV- fc , the ratio products 

A(f) N U*>) N U*i) :_ CS (7) 

are computed and fitted with the functions po ■ (1 + Aj ■ sin $j) to extract the raw Collins 
and Sivers asymmetries. Due to the smallness of the asymmetries involved, Aj is to a good 
approximation Aj = + e ^d + e j,u + e j,d- The fit is performed in the interval (0, 2n) which 
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Figure 5: Difference of the invariant mass of the hadron pair after cuts to the Kg mass. 
The shaded region shows the accepted Kg. 



is divided into 16 bins for charged hadrons and 8 bins for Kg due to the low statistics. 
From Monte Carlo simulations it is found that the angular resolution is much better than 
the bin size. 

The Collins and Sivers asymmetries have been evaluated separately in each kine- 
matic bin and for each data taking period; the results for the different periods have been 
combined using the weighted mean. Finally using the purities a purity matrix P can be 
written 



P 



P. 
P 



here P n ^ (Pr,k) is the fraction of real pions (kaons) inside the identified pion (kaon) 
sample and P n ,K {Pr,-k) is the fraction of misidentifled kaons (pions) inside the pion 
(kaon) sample. Due to the fact that the contribution of other particles (like protons) was 
found to be negligible it is assumed that Pk,k = 1 — P-k,k and Pk,k = 1 — Pr,k- The 
relation between measured A m and corrected A c asymmetries is 



This gives for the kaons 



A c 



p~l A r 



with A 




(9) 



1 



Pn,n + Pr,K — 1 - 



P A m 



A-Pk,k)AI 



(10) 



This equation contains two terms; the first accounts for the dilution of the kaon asymmetry 
due to the purity of the sample, the second removes the contribution of the pion asymmetry 
from the identified kaons. Since both pion and kaon asymmetries on the deuteron target 
are compatible with zero over the full range, this second term has been put to zero, 
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Figure 6: Collins asymmetry against x, z and f° r the "all" charged pions and kaons 
samples from the 2003-2004 data, and the "all" K s 's sample from the 2002-2004 data. 



both in Eq. (flOj) and in the equation giving and the errors have been evaluated 
correspondingly 

The same approach was followed for the neutral kaons, where it has been checked 
that the asymmetry of the background under the Kg peak was compatible with zero. 

6 Results and conclusion 

The final results for the Collins and Sivers asymmetries Acoii and As% v for charged 
pions and charged and neutral kaons on the deuteron target vs. the three kinematic 
variables x, z and j)\ are given in Figs. In the figures, the data points for negative 

hadrons, which are calculated in the same x-, z- and pr-bin as for the positive hadrons, 
have been slightly shifted for graphical reasons. 

Extensive studies to evaluate the size of the systematic error have been performed. 
For some of these studies the z cut has been opened and the data sample has been 
enlarged by a factor of three. The measured Collins and Sivers asymmetries were checked 
against stability among the five different periods of data taking, against the use of different 
estimators to extract the asymmetries, against the reduction of the fiducial volume of 
our spectrometer and against the influence of the trigger system of the experiment. In 
these studies no deviations from the real asymmetries beyond the expected statistical 
fluctuations was observed. Furthermore experimental false asymmetries have been studied 
by combining the data set in such a way that the extracted asymmetries are expected 
to be zero. During all these tests no asymmetries deviating from zero with statistical 

x ) All the numerical values, including the purities, are available on HEPDATA. 
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Figure 7: Sivers asymmetry against x, z and p T for the "all" charged pions and kaons 
samples from the 2003-2004 data, and "all" K%'s sample from the 2002-2004 data. 



significance was observed. 

Also, the correlation between the measured Collins and Sivers asymmetries which 
originates from the non-uniform <fih/<fis acceptance of the spectrometer has been studied 
and the corresponding systematic error has been evaluated to be negligible as compared 
with the statistical error. The smallness of the asymmetries makes the systematic error 
due to the uncertainties on Px and / totally negligible. These studies altogether lead to 
the final conclusion that the systematic errors are considerably smaller (well below 30%) 
than the statistical errors. 

All the measured asymmetries are small, a trend which was already observed in the 
published data of the non-identified hadrons. Small asymmetries are not a surprise, it was 
expected that transverse spin effects be small in the deuteron due to the opposite sign 
which was predicted for the u- and d-quark distributions, very much like in the helicity 
case. 

The interpretation of the results on the deuteron can be done only in conjunction 
with corresponding proton data, measured by the HERMES Collaboration albeit at lower 
energy. Proton target data have been collected by COMPASS in 2007, but the results are 
not final at the time of writing. As shown in Refs. [8,11] a simple analysis of the HERMES 
charged pion data and of the non-identified charged hadron data in COMPASS, assuming 
that all the hadrons are pions, led to the following conclusions: 

1. the favoured and unfavoured Collins functions have about the same size and the 
COMPASS deuteron data are needed for the extraction of the d-quark transversity; 

2. the null result for the Sivers asymmetry for the COMPASS data is a clear indication 
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Figure 8: Collins asymmetry against x, z and p\ for the "leading" charged pions and 
kaons samples from the 2003-2004 data, and "leading" Kg's sample from the 2002-2004 
data. 



that the u- and d-quark Sivers distribution functions have about the same size and 
opposite sign. 

The same conclusions have been obtained in several analyses, using more sophisticated 
tools (see e.g. Refs. [14-16]). A first global analysis which combined the 2002-2004 HER- 
MES pion Collins asymmetries, the COMPASS results for non- identified hadrons, and the 
BELLE data has recently allowed to extract the Collins functions and, for the first time, 
the transversity distributions for the u- and d-quark [17]. Similar analyses can be now 
done including the present pion data which put more stringent constraints. 

The kaon data again show small asymmetries. In the case of the charged kaons, 
although the statistics is a factor of about 6 smaller than for the pions, the error bars are 
still rather small. The neutral kaon sample is smaller in size by a factor of about 3 with 
respect to the charged kaons, and the error bars start being large. The COMPASS data do 
not exhibit the large difference between K + and 7r + asymmetries seen by HERMES. Very 
much like for the 7r =l= case, cancellations are expected between u- and d-quarks when using 
the isoscalar deuteron target. Therefore the smallness of the COMPASS kaon asymmetries 
suggests that the sea quark contributions to the asymmetries are small. The kaon data 
provide a unique handle on the s-quark, but in this case the sea-quark contributions can 
not be neglected, and a full global analysis including pions and kaons is mandatory. 

To summarise, COMPASS has made the first precise measurements of the Collins 
and Sivers asymmetries for charged pions with a transversely polarised deuteron target. 
The same asymmetries have also been obtained for charged and neutral kaons. All the 
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Figure 9: Sivers asymmetry against x, z and ~p\ for the "leading" charged pions and kaons 
samples from the 2003-2004 data, and "leading" K° s 's sample from the 2002-2004 data. 



measured asymmetries are small, pointing at a cancellation between the u- and d-quarks 
contributions. More quantitative information, in particular for the s-quark distributions, 
can be obtained with global analyses, in which the COMPASS measurements with a 
transversely polarised proton target undoubtedly will play an important role. 

We acknowledge the support of the CERN management and staff, the special sup- 
port of CEA/Saclay in the target magnet project, as well as the skills and efforts of the 
technicians of the collaborating institutes. 
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